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Abstract 

We investigate the feasibility of using laser accelerated protons/deuterons 
for positron emission tomography (PET) isotope production by means of the 
nuclear reactions ^^B(p, n)^^C and ^°B((i, n)^^C. The second reaction has a 
positive Q-value and no energy threshold. One can, therefore, make use of 
the lower energy part of the laser-generated deuterons, which includes the 
majority of the accelerated deuterons. The ^^C produced from the reaction 
^''B((i, n)^^C is estimated to be 7.4 x 10^ per laser-shot at the Titan laser at 
Lawrence Livermore National Laboratory. Meanwhile a high-repetition table 
top laser irradiation is estimated to generate 3.5 x 10'' ^^C per shot from the 
same reaction. In terms of the ^^C activity, it is about 2 x 10^ Bq per shot. 
If this laser delivers kHz, the activity is integrated to 1 GBq after 3 minutes. 
The number is sufficient for the practical application in medical imaging for 
PET. 



1. Introduction 

Laser accelerated proton beams have been applied to produce positron 
emitting isotopes by several groups [H El 13 S]. Those groups investigate 
the feasibility of ^^C isotope production rate by the laser-induced reaction 
^^B(p, n)^^C. The reaction has a negative Q-value of -2.765 MeV and, there- 
fore, requires large scale lasers like the VULCAN laser facility at the Ruther- 
ford Appleton Laboratory (RAL). Such large scale lasers are single shot. At 
current technology table-top lasers, which generate low temperature plasma 
but operate at high repetition rate, cannot lead to the reaction. In Ref. [1] 
from the reaction ^^B(p, n)^^C induced by the laser-generated protons and 
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the secondary boron target, the activities of 2 x 10^ Bq of the isotope ^^C 
have been measured. They have carried out the experiments using the VUL- 
CAN laser, which dehvers pulses of energies up to 120 J, duration 0.9-1.2 ps 
and a wavelength of 1053 nm. This laser pulse focused on a thin Al foil gen- 
erates protons of energies as high as 37 MeV, from the contaminants on the 
foil target. Lately, the experiment is repeated at the same VULCAN laser 
at a pulse energy of 300 J with a pulse duration of 750 fs [3], a ^^C activity 
of 6 X 10^ Bq has been reported. On the other hand the activity required 
for PET is at least 200 MBq and ideally it is 800 MBq p. 

In this paper we investigate the feasibility of this application but us- 
ing the reaction ^^B{d,ny^C, instead. This reaction has positive Q-value. 
There is, therefore, no energy-threshold for the reaction |5]. This means that 
smaller size lasers with higher-repetition rates can be used in contrast to the 
large-scale lasers, though such lasers generate proton beams at lower energies. 
Laser generated deuterons from deuterated plastic layer on an Al foil target 
has been investigated only recently [HI El El E]- A numerical investigation 
suggests that by eliminating surface contaminants, which contain protons, 
deuterons can be accelerated successfully, even though the laser pulse energy 
is relatively low. We compare estimates of ^^C yields for the proton- and 
deuteron-induced reactions irradiated by two types of lasers. First the Titan 
laser at the Jupiter Laser Facility at Lawrence Livermore National Labora- 
tory (LLNL) [TD]. The laser pulse delivered by this laser has energies up to 
150 J, the pulse length of 0.6 ps and the wave length 1054 nm ^U\. These 
are close to the parameters of the VULCAN laser at the configuration in 
Ref . [1] . The repetition rate of the Titan laser is 2-shots per hour. The other 
is a table-top laser, ASTRA Ti:sapphire laser facility [TT] at RAL, with a 
repetition rate of 2 Hz, a laser energy of 200 mJ, a pulse length of 60 fs and 
a wave length of 790 nm. 

2. Ion acceleration by lasers 

Several mechanisms are proposed to be responsible for the ion accelera- 
tion by interaction of high-intensity laser pulses with thin foil targets. The 
target normal sheath acceleration (TNSA) is one of those mechanisms and 
is thought to be the dominant mechanism at laser intensities 10^^ to 10^^ W 
cm~^. We assume the TNSA mechanism as the mechanism of the proton 
and deuteron acceleration. In the electric field of a laser, electrons on the 
front surface of the target quiver. The oscillating motion of the electrons 
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causes a force, which drives the electrons into the target normal direction, 
due to the oscillational magnetic field perpendicular to the electrons' motion. 
The accelerated electrons penetrate through the target, ionizing the surface 
contaminants on the back, and form a dense sheath on the back. The charge 
separation between the escaping electrons and the ions left behind causes 
an electrostatic field which accelerates the protons in the hydrocarbon con- 
taminant on the back surface. Under the TNSA mechanism, the protons 
and deuterons are expected to be accelerated by the same force eE, where 
e and E are the electron charge magnitude and the electric field due to the 
charge separation, respectively. The energy of the deuterons accelerated by 
the electrostatic field is determined by the charge of the deuteron which is 
the same of protons. Consequently the energy of the deuterons should be the 
same as the energy of the protons accelerated by the same field, while their 
velocities are different. On this account we expect that the deuteron energy 
spectra is similar to that of the protons. 

However, recently published experimental data [SI El E] show a controver- 
sial result. Experiments on the deuteron acceleration have been performed 
using T-cubed laser at the university of Michigan. This laser delivers a pulse 
energy up to 6 J in a 400 fs. The target was a 6 fim thick Al or 13 /zm 
thick Mylar, coated with a layer of deuterated polystyrene either on the back 
or on the front or both side of the target. For the detection of accelerated 
ions, a Thomson parabola spectrometer has been used with a micro channel 
plate with a phosphor screen. The experimental observation that protons 
are accelerated preferentially to deuterons [HI |8] is because the protons in 
the surface contaminants are on the top of the deuterated plastic layer. The 
protons, which are facing directly to the charge-separation electric field, see 
the maximum field strength, while the deuterons underneath the surface 
contaminant see the reduced electric field shielded by the protons. If the 
electrostatic field formed on the back surface is large enough, the layer of 
protons are ablated totally by the field. Then the ions of the target become 
exposed to the field and are accelerated by the remaining electrostatic field 
energy. Another reason could be that if you have a mixture of different par- 
ticles, the ones with lower mass leave the target first because of the higher 
velocity. Thus once they leave the target the electric field is reduced and 
the heavier ones see a lower field. Thus it is important to have pure targets 
to maximize the deuteron fiux. The surface contaminants layer can be re- 
moved either by thermal heating [H] or using a secondary laser to ablate the 
surface. The use of high-repetition lasers has an advantage from this view 
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point, since the contaminants of lower mass (protons) will be removed in the 
repetitional pulses. The removal of the contaminant allows acceleration of 
the target material (deuterons) efficiently [71 19] . 

The reductions both in the maximum deuteron energy and in the num- 
ber of accelerated deuterons compared to the case of protons is observed 
to be approximately 1:10 [S] at the given experimental condition, i.e., the 
thicknesses of the target and of the contaminant layer, and the specific laser 
parameters. Further experimental investigations are needed to understand 
difference between the proton- and deuteron-spectra, generated at different 
laser parameters and experimental conditions. For the purpose of the yield 
estimate of positron emitting nuclei, we use the same energy spectra as pro- 
tons for deuterons at a set of given laser parameters. We stress that this 
alternative reaction to produce radioisotopes using deuterons is in any case 
interesting, considering the possibilities that with advancing technologies, 
higher repetition lasers of relatively low energy and short pulse might be 
feasible. 



3. Yield and angular distribution of neutrons 

3.1. the reaction ^^B(p,n)^^ C 

The yield from the reaction ^^3(^,72)^^0 is estimated at Titan laser ir- 
radiation at LLNL. By irradiating Au foil target of thicknesses from 5 fim 
to 250 fim, coated with 200 nm of ErHa on the rear surface, proton plasma 
which has the highest energy 36 MeV have been observed [10]. The proton 
spectra from this laser irradiation interacting with a 14 nm thick Au target 
coated with ErHs is measured with RC film. The spectra are parametrized 
by the following two temperatures distribution: 

where Ti = 3.3 MeV, Ta = 13.5 MeV, Ni = 1.2 x lO^^ ^nd N2 = 2.3 x 
10^^ [T2I Ho]. We reproduce the observed proton spectra in the top panel in 
Fig. [Tj If one locates a secondary boron slab target after the Au target, ^^B 
in the target undergoes the nuclear reaction ^^B(p, n)^^C. In this simulation, 
we assume that the secondary target is a 3 mm-thick slab of natural boron. 
The thickness of the target is chosen to be same as in the experiment at 
VULCAN laser [1]. We call this target (a). Given that the proton spectra 
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spread as high as 36 MeV, the secondary target must be thick to ensure 
that all the accelerated protons which go through the boron target cause the 
nuclear reaction. 

The ^^C yield (Nf) from the laser irradiation is evaluated by [5] 

Nf = Nia{E)nTd, (2) 

where Ni,a{E),nT and d are the number of accelerated ions, the cross sec- 
tion of the reaction, the density and the effective thickness of the boron 
target, respectively. The cross section for the reaction ^^B{p,ny^C is shown 



in Appendix in Fig. |A.4| by triangles and is roughly 100 mb at the energy 
higher than 5 MeV. The target density is nT=8 x 10^^ cm~^. We replace 
the effective thickness of the target by the projected range of the protons in 
the target |13]. The projected range is given in a table as a function of the 
proton incident energy. For a rough estimate of the yield, we approximate 
the proton energy in the plasma, which spread from to 36 MeV in reality, 
by the most probable energy [H] for the reaction. The most probable energy 
for the plasma ions that cause the nuclear reaction is given by the Gamow 
peak energy (Eq) [I5] and is found using the saddle point method, i.e.: 

d / E 



, 1 , 



where b = 3I.28Z1Z2M2 (keVa ), denoting the atomic number of the colliding 
nuclei Zi, Z2 and the reduced mass number M. The reduced mass number 
M = AiA2/{Ai + A2), where Ai and A2 are the mass numbers of the colliding 
ions. The temperature kT is replaced by the plasma temperature, then. 




Eg=^] . (4) 



For the case of proton plasma generated by the Titan laser Eq is 3.9 MeV. 
At this energy the cross section of the reaction and the projected range of 
the protons are about 50 mb and d =150 yum, respectively. A rough estimate 
of the fusion yield for the Titan laser facility gives 

Nf = 10^^50(m6)8 x 10^^{cm-^)150{fim) ~ 5. x 10^ (5) 

per shot. We note that the effective energy used in the estimate is higher 
than the threshold of the reaction. 
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For a more precise estimate of the yield, we calculate the angular distri- 
bution of the reaction products, i.e., neutrons, from the given proton spectra. 
From the kinematical relations the emission angle of ^^C is limited to less 
than 60 deg. The energy of generated ^^C is also limited to less than 3 MeV. 
Given that the emission angle is forward and the energy is very low, all ^^C 
are stored in the target. The thickness of the secondary target is enough to 
let ^^C lose its energy inside the target. Experimentally the activity of ^^C is 
measured by detecting 0.51 MeV gamma-rays by Nal scintillation detectors. 
For completeness we determine ^^C yield from the neutron yield, which could 
bring other informations on the reaction mechanism in the plasma. 

The angular distribution of the neutron yield is assessed by 

'''' -m,mnr4E), (6) 



OEdn ' dEdVL 

where </>(-£') and are the proton energy spectra and the differential 

cross section of the reaction; the proton spectra is parametrized in Eq.Q 
and shown in the top panel in Fig. [l] d{E) is the projected range which is a 
function of the incident energy, as mentioned before. We use the assumption 
of an isotropic angular distribution of neutrons in the center-of-mass (CM) 
system. The differential cross sections at different angles are obtained by 
transforming the cross sections from the CM system to the laboratory system. 
In Fig. [T] the lower panels show the neutron spectra at emission angles, 0, 
45, 90, 135 and 175 degrees, from the proton energy distribution in the 
top panel. One can see that neutrons have broad energy spectra from to 
30 MeV, especially at the forward angle. 

By integrating the neutron yield, the total yield of the reaction ^^B(p, n)^^C 
is estimated to be 1.5 x 10^ at Titan. This is in good agreement with the rough 
estimate given above. 

3.2. The reaction ^^B(d,n)^^C 

^^C ions are generated by the deuteron-induced reaction ^'^B((i, n)^^C, as 
well. This reaction has a positive Q-value 6.5 MeV. If one uses deuteron 
enriched hydrocarbon coating on a Au foil target and replace the secondary 
natural boron target by a ^'^B enriched target, ^°B undergoes the nuclear 
reaction ^°B((i, ?t,)^^C. Experimental data of deuteron spectra at such a setup 
is not yet available. Assuming that the deuteron energy spectra is the same 
as the proton energy spectra, we estimate the yield for ^°B((i, n)^^C similarly 
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Figure 1: Proton spectrum at Titan laser irradiation on a An target coated with ErHs 
(top panel). Expected neutron spectra for the reaction ^^B(p, n)^^C on a 3 mm-thick 
target at emission angles, 0, 45, 90, 135 and 175 degrees, at Titan laser irradiation. 
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to the reaction ^^B(p, ?2)^^C. Eq. ^ gives again a rough estimate of the 
fusion yield by substituting A^j=10^^, a = 200 mb, ut = 8 x 10^^ cm~^ 
and d = 150 /im, one obtains Nj ~ 2 x 10^. In Eq. (|6| by substituting the 
differential cross section by that of the reaction ^^B{d,n)^^C, we estimate 
the neutron yield from the interaction between the deuteron plasma and the 
secondary enriched target. Fig. [2] shows the neutron yield estimated for 
the reaction ^°B((i, n)^^C at 5 selected emission angles. The neutron spectra 
range from 5.93 MeV, due to the reaction Q-value, to more than 20 MeV at 
all the detection angles. We note that the secondary boron target is chosen 
to be 5 mm thick to ensure that all the deuterons which go through the 
boron target cause the nuclear reaction. By integrating the neutron yield, 
the total yield of the reaction ^^B{d,ny^C is estimated to be 7.4 x 10^ 
per shot at Titan. This is again in good agreement with the former rough 
estimate. We also notice that neutrons could be produced by the deuteron 
disintegration. We expect this effect to be small because of the Q-value, but 
it could be precisely determined from kinematical considerations and from a 
precise measurement of the neutron yield and the gamma decay from ^^C. 

Spencer et al. [H] studied proton spectra from plastic, aluminum and 
copper targets from a high-repetition-rate tabletop-laser irradiation. The 
laser parameters in the experiment are the laser energy of 200 mJ, the pulse 
length of 60 fs and the wave length of 790 nm. It is of interest that the 
plastic target irradiation, at such laser parameters, generates higher energy 
protons and that the proton spectra have characteristic shape which cannot 
be approximated by a simple Maxwellian [T3]. Thus we choose a proton 
spectra from Mylar target irradiation, which are shown in Fig. 2 in Spencer 
et al. [11] , and fit the experimental data by the following formula: 



<P{E) = c—^ ^ exp{-E/kTHH) 
where E is the energy of the protons in the laboratory system and 



E' = E- 2^EEo + Eo; (8) 

c,Ci,Thh,Ti and Eq are fitting parameters; k is the Boltzmann constant. 
Eq represents the energy of an extra acceleration as discussed in a previous 
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Figure 2: Expected neutron spectra from the deuteron-induced reaction ^''B(d, n)^^C on 
the secondary 5 mm-thick ^"B enriched target at 5 selected emission angles at Titan. 



9 



work [H]. For the spectra of the irradiation on 36 /xm-thick target, which has 
the number of accelerated ions larger than the other spectra, the obtained 
fitting parameters are listed in Tab. [1} The proton spectrum given by these 
parameters is reproduced in the top panel in Fig. [3j Assuming that the 
deuteron spectra in deuterated plastic target irradiation are similar to the 
proton spectra, we estimate the total yield for the reaction ^^B{d, nY^C. The 
angular distribution of the neutron yield with the secondary boron target is 
calculated by Eq. ([6]). The neutron spectra at 5 selected emission angles are 
estimated and shown in Fig. |3| The neutron spectra range from 5.9 MeV 
to 8.1 MeV at the forward angle and to 6.4 MeV at the backward angle. 
The evident two peaks in the neutron spectra are due to the two peaks in 
the proton spectra in the top panel. By integrating the neutron yield, the 
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Table 1: Fitting parameters for 36 /im Mylar target irradiation in Spencer et al 



total yield of the reaction ^^B{d, nY^C is estimated to be 3.5 x 10^ per shot. 
While in the reaction ^^B(p, n)^^C there will be no ^^C yield at this table-top 
high-repetition laser irradiation. 

4. Discussion 

The total yield of the reaction ^^B{p,ny^C is estimated to be 1.5 x 10^ 
per shot at Titan laser irradiation. While the total yield of the reaction 
^^B[d,nY^C is estimated to be 7.4 x 10^ per shot at the same laser irradi- 
ation. The estimated yield of the reaction ^^B{d,ny^C is five times higher 
than that of the other reaction, given the same laser-pulse irradiation condi- 
tion. This yield gain by changing the reaction is explained by the fact that 
the deuteron-induced reaction has a positive Q- value, while the other proton- 
induced reaction has a negative Q-value. That is a large part of the plasma 
ions, which have energy less than 3 MeV, do not contribute to the yield in 
the reaction ^^B{p,ny^C. If we use a high repetition-rate tabletop-laser, the 
yield is estimated to be 3.5 x 10''' per shot from the reaction ^^B{d, n)^^C. We 
note that the yield for the reaction ^^B{d,nY^C could be overestimated due 
to our assumption of the identical energy spectra for protons and deuterons. 
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Figure 3: Proton spectrum at a table-top laser irradiation on a 36 /xm-thick Mylar target 
(top panel). Expected neutron (red curves) and ^^C (blue curves) spectra from the reaction 
^'^'B{d,n)^^C on 3 mm-thick enriched target at 5 selected emission angles at a high 
repetition-rate table-top laser irradiation. 
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If the major part of ions accelerated by this laser have energies less than 
0.5 MeV, while the cross section of the reaction ^^B{d,ny^C is small at 



the incident energy less than 0.5 MeV, as shown by diamonds in Fig. A. 4 
the amount of ^^C yield might be increased by using a more intense high- 
repetition laser. Here we assume a high-repetition laser, which meets the 
above requirement, at a laser energy 1 J and a pulse duration 40 fs and a 
repetition rate 10 Hz, and estimate the yield from this laser. These laser 
parameters are, indeed, similar to the parameters of the high repetition laser 
in Fritzler et al. [2] . By studying laser-driven proton scaling laws [161 113 IH] , 
we estimate the plasma temperature of 10 MeV and the number of accelerated 
ions of 10^^ for this laser. Moreover we assume the maximum proton energy 
of 10 MeV. For the proton spectra characterized by these parameters, the 
yield for the reaction ^^B(p, n)^^C and for ^°B((i, 72)^^0 are estimated to be 
1.3 X 10^ and 2.8 x 10^, respectively, per shot. 

Spencer et al. [Ij and Fritzler et al. [2] discuss the application of the 
laser generated proton plasma to the PET isotope production. Spencer et 
al. [1] used the reaction ^^B(p, n)^^C at VULCAN laser, while Fritzler et 
al. [2] used a table top laser, at a laser energy of 840 mJ and a pulse length 
of 40 fs, and at a repetition rate of 10 Hz. Spencer et al. observed ^^C 
activities of 2x 10^ Bq per shot and estimates 10^ Bq activities after 500 s(8.3 
minutes) of irradiation, assuming 10 Hz operation of VULCAN laser. Fritzler 
et al. [2] estimated the activity 13.4 MBq after 30 min irradiation at the laser 
at Laboratoire d'Optique Appliquee (LOA). 

The result in Spencer et al. [T] could be compared with our estimate 
at Titan laser irradiation (1.5xl0^x 0.69/20.4/60=850 kBq) using the same 
reaction. The difference comes from the fact that proton energy spectra 
from Titan laser-irradiation (Top panel in Fig. [T| has a higher flux than 
that from VULCAN laser-irradiation (Fig. 4 in Ref. [Ij). The increase of the 
proton flux could be caused by the usage of the Au target coated with ErHs 
in the experiment at LLNL, while protons originate only from impurities 
on the target in the experiment at VULCAN. In the case of the reaction 
^°B((i, n)^^C at Titan laser irradiation the produced ^^C activity is estimated 
to be (7.4xl09x 0.69/20.4/60=4200 kBq). The estimated activity is still 
much smaller than the required activity for PET isotope, i.e. 800 MBq. 

If one uses the reaction ^^B{d,ny^C and the table-top, high-repetition 
laser, we have estimated ^^C production of 3.5 x 10^ per shot, in terms of the 
"C activity it is (3 xlO^ x 0.69/20.4/60 ~ 2 x 10^ Bq) per shot. Considering 
that ^^C decays with the half-life of 20.4 min and given 2 Hz repetition-rate. 
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after 30 minutes of irradiation the activity is integrated to be 22 MBq. If 
this laser dehvers kHz, the activities are integrated to 1 GBq after 3 min- 
utes and meet the required activity for the practical use of the PET isotope 
production. 

Alternatively, for another high-repetition laser which delivers higher en- 
ergy pulses of 1 J and a pulse duration of 40 fs and a repetition rate of 
10 Hz, the ^^C activities for the reactions ^^B(p, n)^^C and ^°B((i, n)^^C 
are estimated to be (1.3 xlO^ x 0.69/20.4/60 ~ 730 Bq) and (2.8 xlO^ x 
0.69/20.4/60 ~ 1.5 kBq), respectively, per shot. Given 10 Hz repetition-rate, 
after 30 minutes of irradiation the activities are integrated to be 8.2 MBq 
and 17 MBq, respectively, for the reactions ^^B(p, n)^^C and ^°B((i, ?2)^^C. 

We remind that the ^^C activities from the reaction ^'^B((i, nY^C at given 
laser parameters are estimated under the assumption of identical proton- 
and deuteron-energy spectra. Taking into account the recent experimental 
observation |8], the ^^C activities from the reaction ^^B{d,nY^C are likely to 
be reduced than the above estimates. 

5. Summary 

We have studied the feasibility of ^^C isotope production by means of the 
nuclear reactions ^^B(p, n)^^C and ^^B{d,nY^C induced by laser-accelerated 
protons/deuterons. We have estimated ^^C produced by the reaction ^^B(j9, riY^C 
to be 1.5 X 10^ per shot at Titan laser irradiation. While ^^C produced by 
the reaction ^°B((i, 71)^^0 at the same laser irradiation is estimated to be 7.4 
X 10^. The latter deuteron-induced reaction increases the yield 5 times than 
that of the proton-induced reaction. However the yield per shot is still too 
low for the practical use in the PET isotope production. 

The yield of ^^C at a table top, high-repetition laser has been estimated 
to be 3.5 X 10^ per shot for the deuteron-induced reaction. In terms of the 
^^C activity, it is about 2 x 10^ Bq per shot. If this laser delivers kHz, the 
activity is integrated to 1 GBq after 3 minutes and meet the required activity 
for the practical use of the PET isotope production. 

The ^^C activities for the reaction ^^B{d,ny^C are likely to be overes- 
timated due to the assumption of the identical energy spectra for protons 
and deuterons. The realization of this potential application requires further 
research to understand optimum parameters for deuteron acceleration by in- 
tense lasers. Our results have demonstrated a promising alternative way of 
using high-intensity laser for the application in the PET isotope production. 
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Figure A. 4: Experimental data of cross sections for the reactions ^^B{p, ?t,)-'^-'^C [TS] (trian- 
gles) and ^°B((i, n)^^C [20j (diamonds) are shown as functions of the incident energy both 
in the center-of-mass system. 
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Appendix A. Cross sections 

Experimental data of cross sections for the reactions ^^B(p, 'ri)^^C and 
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B((i, n)^^C are compared in Fig. A. 4 
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